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Abstract 

In three dimensions, the gluon condensate of pure SU(3) gauge theory has ultraviolet diver- 
gences up to 4-loop level only. By subtracting the corresponding terms from lattice mea- 
surements of the plaquette expectation value and extrapolating to the continuum limit, we 
extract the finite part of the gluon condensate in lattice regularization. Through a change of 
regularization scheme to MS and (inverse) dimensional reduction, this result would determine 
the first non-perturbative coefficient in the weak-coupling expansion of hot QCD pressure. 
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1. Introduction 



As non-Abelian gauge theories in three and four dimensions are confining, their proper- 
ties need, in general, to be studied non-perturbatively. If the observables considered in- 
volve momenta or masses (M) large compared with the confinement scale, however, then the 
conceptual framework of the operator product expansion may allow to isolate the non- 
perturbative dynamics into only a few (gluon) condensates, while the rest of the answer can be 
computed by perturbative means. A classic example is the mass of a heavy quark-anti-quark 
bound state P|. The task faced by numerical lattice simulations might then get significantly 
simplified, as local condensates are simpler to measure with controlled systematic errors than 
correlation functions of heavy states. 

On the other hand, the physical interpretation of a "bare" lattice measurement of a gluon 
condensate is non-trivial. The reason is that the condensate is represented by the expectation 
value of a dimensionful singlet operator and, in general, contains ultraviolet divergences of the 
same degree as its dimension. Operator product expansion type relations are often derived 
employing dimensional regularization, since the system then only contains one large param- 
eter (M) rather than two (M and the momentum cutoff), which simplifies the derivation 
considerably. Making use of lattice results in such a context requires then a transformation 
from lattice to continuum regularization. While in principle a well-defined perturbative prob- 
lem (see, e.g., Refs. this is in practice somewhat problematic in four dimensions, given 
that there are contributions from all orders in the loop expansion. 

The observable we consider in this paper is the (lowest-dimensional) singlet gluon conden- 
sate in three dimensions (3d), measured with pure SU(3) gauge theory. As 3d pure Yang- Mills 
theory is super-renormalisable, the problem of changing the regularization scheme becomes 
solvable: there are ultraviolet divergences up to 4-loop level only 4 . Furthermore, as we will 
elaborate in the following, all the divergences have been determined recently. These advances 
allow us to obtain a finite "subtracted" continuum value for the gluon condensate in lattice 
regularization. A conversion to the MS scheme, amounting to the (perturbative) computation 
of the constant C4 in Eq. ()2.12|) below, remains however a future challenge. 

There might be various physics settings where the 3d gluon condensate finds practical 
applications. The one that motivated us, is that this condensate appears in high-temperature 
physics, where the temperature T plays the role of the mass scale M mentioned above. Indeed 
3d pure Yang-Mills theory determines the leading non-perturbative contribution to the weak- 
coupling expansion of the pressure (and a number of other quantities) of physical QCD 012], 
through a conceptual counterpart of the operator product expansion, called finite-temperature 
dimensional reduction |3IHlini- Other applications might exist as well. 

The plan of this paper is the following. In Sec. |21 we specify the observables considered 
and discuss the theoretical setting of our study. Numerical results from lattice Monte Carlo 
simulations are reported in Sec. 121 and we conclude in Sec. ^ 
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2. Theoretical setting 



We start this Section by formulating the observables that we are interested in, in the formal 
continuum limit of the theory. The ultraviolet (UV) divergences appearing in loop contribu- 
tions are at this stage regulated through the use of dimensional regularization. Later on we 
go over to lattice regularization, in order to give a precise non-perturbative meaning to the 
observables introduced, allowing for their numerical determination. 

The Euclidean continuum action of pure SU(A'c) Yang- Mills theory can be written as 



Se 



Jd^xCE, CE = ^Jr[Fii]. (2.1) 



Here d = 3 - 2e, is the gauge coupling, k,l = I, ...,d, Fki = i[Dk,Di], = dt - iA^, 
Ak = AIT", are the Hermitean generators of SU(iVc), normalised as Tr [TT''] = (5"V2, 
and repeated indices are assumed to be summed over. Leaving out for brevity gauge fixing 
and Faddeev-Popov terms, the "vacuum energy density" reads 



/ms = - T-r In 



j VAk exp(^-5E 



(2.2) 



where V is the d-dimensional volume, T>Ak a suitable (gauge-invariant) functional integration 
measure, and we have assumed the use of the MS dimensional regularization scheme to remove 
any 1/e poles from the expression. We note that fj^ has the dimensionality [GeV]*^. 

In strict dimensional regularization, f-^ of course vanishes order by order in the loop expan- 
sion, due to the absence of any mass scales in the propagators. This behaviour is unphysical, 
however, and due to an exact cancellation between UV and infrared (IR) divergences; for 
an explicit discussion at 3- loop level in a related case, see Appendix B of Ref. [I^. In fact 
non-perturbatively the structure of is rather 



/ms — 53 



(4vr) 



12 768 ) 2N^gl 



(2.3) 



where dA = — 1, and we have introduced an MS scheme scale parameter /x. The coefficient 
of the logarithm in Eq. (|2.3|) has been determined in Ref. with a perturbative 4-loop 
computation, by regulating all the propagators by a small mass scale niQ, and sending niQ 
only after the computation (see also Ref. JU). The non-perturbative constant part Bq,^ 
which actually is a function of Nc, is what we would ultimately like to determine. 

One direct physical application of Bq is that it determines the first non-perturbative contri- 
bution to the weak-coupling expansion of the pressure p of QCD at high temperatures [21 E]- 
To be precise, this contribution is of the form 6p = dAN^ g^T'^ Bq / (At:)^ , where g"^ is the renor- 
malised QCD gauge coupling. Terms up to order 0{g^ ln{l/g)) are, in contrast, perturbative, 
and all known by now jl2j . 

^In Ref. |12| . Bq was denoted by /3q, but we prefer to introduce a new notation here, in order to avoid 
confusion with the coupling constant l3 appearing in Eq. H2.6|l . The subscript G might refer to gluons. 
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For future reference, we note that given /j^, we immediately obtain the gluon condensate: 

(2.4) 



1 

25i 



2 d 

-93 



35! 



■,dANl 
(47r)4 



43 
12 



157 
768' 



In. 



^ i)+i?G + 0(e) 



(2.5) 



We now go to the lattice. In the standard Wilson discretization, the lattice action, Sa 
corresponding to Eq. (|2.H) . reads 



Sa 



1 



X k<l 



(2.6) 



where Pfcz(x) = C/fc(x)C//(x + fc)?7^ ^(x + /)C/j ^(x) is the plaquette, C/fc(x) is a link matrix, 
X + = X + aik, where a is the lattice spacing and ik is a unit vector, and 



2A^c 



(2.7) 



Note that the gauge coupling does not get renormalised in 3d, and the parameters ap- 
pearing in Eqs. (|2.H) . (|2.7j) can hence be assumed finite and equivalent. The observable we 
consider is still the vacuum energy density, Eq. (|2.2|) . which in lattice regularization reads 



fa 



lim — In 

y^oo V 



jvUk exp(^-5'a 



(2. 



where T>Uk denotes integration over link matrices with the gauge-invariant Haar measure. 

Now, being in principle physical quantities, the values of and fa must agree, provided 
that suitable vacuum counterterms are added to the theory. Due to super-renormalizability, 
there can be such counterterms up to 4-loop level only and correspondingly 



A/ 



fa fms 



(2.9) 



093 



C7i ^( In^ + C( ) + C2 % + C3 ^ + C453'( In^ + + 0{gla) , (2.10) 



93 



where the Ci are dimensionless functions of Nc- The values of Ci,C2,C3,C4 are known, as 
we will recall presently; C{ is related to the precise normalisation of the Haar integration 
measure and void of physical significance; and C4 is unknown as of today. 

Correspondingly, the gluon condensates, i.e. the logarithmic derivatives of f-^, fa with 
respect to 53, can also be related by a perturbative 4-loop computation. Noting that three- 
dimensional rotational and translational symmetries and the reality of Sa allow us to write 



2 d f 



^(1 



■Tr [P, 



12 



(2.11) 
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and employing Eqs. (\2.5\i . ()2.1U|) . we obtain finally the master relation 



(4vr) 



The values of the constants ci, C4 are trivially related to those of Ci, C4 in Eq. (|2.1()j) . 
Now, a straightforward 1-loop computation yields 

ci = ^ f« 2.66666667 , (2.13) 

where the numerical value applies for A'c = 3. 

The 2- loop term is already non-trivial: it was first computed in four dimensions in Ref. |13j . 
and in three dimensions in Ref. jl4j . The 3d result can be written in the form 

2 dAN^ / 47r2 S2 vr^ 2 \ , , , , 

« = - 3 (4# + T - - T + + rO = ■ 

where the coefficients S, ki can be found in Refs. |SE1) aiid we have defined 

1 /■'^/^ o o y],- sin2j;,sin2(x, •+ yjOsin^u,- , , , , 

«5 = -4 / d^xd^y ^ . , ^- 2 = 1-013041(1) . (2.15) 

The numbers in parentheses in Eqs. 1)2. 14(1 . (|2.15p indicate the uncertainties of the last digits. 

The 3- loop term is well known in four dimensions since a long time ago |16| , but the same 
computation has been carried out in three dimensions only very recently |17j : 

C3 = 6.8612(2) . (2.16) 

This value improves on a previous estimate C3 = 6.90io!i2 ^M, obtained through the evaluation 
of the 3- loop graphs with the method of stochastic perturbation theory jl9j . 

The value of C4 follows by a comparison of Eqs. (|2.3j) and (|2.1fl)) : there is no /i-dependence 
in fa, so that the one in f-^ determines the coefficient of the logarithm in A/. Consequently, 

(IaN^ /43 157 o\ 

The knowledge of ci , C2 , C3 , C4 allows us to subtract all the divergent contributions from 
the gluon condensate. A finite 4- loop term, parametrised by C4 in Eq. (|2.12() . however still 
remains. It could in principle be determined by extending either the method of Ref. or 
of Ref. IIHI to 4-loop level. There is the additional complication, though, that intermediate 
steps of the computation require the use of an IR cutoff, which then cancels once the lattice 
and MS results are subtracted, in Eq. (|2.1fl|) . This computation has not been carried out yet, 
and therefore we will not be able to determine Bq in this paper. We can determine, how- 
ever, the non-perturbative input needed for it (cf. Eq. (|3.1j) below), the purely perturbative 
determination of C4 then remaining a future challenge. 
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p 


volumes 






12 


243, 32^, 48^ 






16 


243, 32^, 48^, 64^ 






20 


243, 32^, 48^ 






24 


{12^, 143, 16^, 20^, 243}, 32^, 


483, 


643 


32 


{143, 16^, 20^, 243, 323}, 433^ 


643, 


963 


40 


{323}, 483^ 643^ 9g3 






50 


{203, 243, 263, 283, 323, 4^3 j.^ 


643, 


963, ;^283, 3203 


64 


{483, 643}, 9g3^ ;l283, 3203 






80 


{643}, ;l283, 3203 






[100 


1283, 1923, 3203] 






[140 


{1283}, 1923, 3203] 






[180 


{1923}, 3203] 







Table 1: The lattice spacings (parametrised by (5, cf. Eq. H2.7|) ') and the volumes (in lattice 
units, A^3, go that V = N^a?) studied. On each lattice we have collected ~ 10^. ..10^ indepen- 
dent measurements. The lattices in curly brackets have been left out from the infinite-volume 
extrapolations, while for the lattices in square brackets the significance loss due to the ul- 
traviolet subtractions in Eq. (|2.12jl is so large (six orders of magnitude or more) that the 
subtracted values have little effect on our final fit (see below). 

3. Lattice measurements 

The goal of the numerical study is to measure the plaquette expectation value, (1 — |Tr [-Pi2])a) 
as a function of /?, such that the extrapolation in Eq. 1)2. 12() can be carried out. For each 
/3, the infinite-volume limit needs to be taken. Given that the theory has a mass gap, we 
expect that finite-volume effects are exponentially small, if the length of the box L is large 
compared with the confinement scale, ~ l/^i- Writing L = Na, where N is the number of 
lattice sites, the requirement L » 1/5(3 converts to /3/N ^ 6 (cf. Eq. (|2.17|) ). Detailed studies 
with other observables show that in practice the finite-volume effects are invisible as soon as 
P/N < 1 |20j . The values of /3 and that we have employed are shown in Tabled Earlier 
lattice measurements of the same observable were carried out with a volume A^3 = 32^, with 
values of /3 up to /3 = 30 

It is important to stress that the subtractions in Eq. 1)2. 12() lead to a major significance 
loss. Essentially, we need to evaluate numerically the fourth derivative with respect to (3^^ 
of the function (1 — ^Tr [Pi2])ai at the point (3~^ = 0. Another way to express the problem is 
that as the numbers ci, ...,C4 are of order unity (cf. Eqs. (|2.13|) ~ H2.17() '). the dominant term, 
C1//3, is about six orders of magnitude larger than the effect we are interested in, ~ 1//3'*, if 
/? ~ 100. Therefore the relative error of our lattice measurements should be smaller than one 
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Figure 1: The plaquette expectation value, "plaq" = (1 — |Tr [^121)0! as a function of 1//3. 
Statistical errors are (much) smaller than the symbol sizes. The dotted curve contains the 
four known terms ci//3 + C2//3'^ + C3 / (3^ + c^ln (] / P'^ from Eq. H2.12() . together with terms of 
the type 1/(3^ and l/f3^ with fitted coefficients. 

part in a million. We also need to know the coefficients Cj with good precision. 

Lattice-measured values of (1 — |Tr [Pi2])a are shown in Fig. ^ as a function of In 
order to demonstrate the accuracy requirements we are faced with, Fig. |2 shows — 
|Tr [i-i2])a, before and after the various subtractions. It is observed from Fig. Qthat after all 
the subtractions, this function indeed appears to have a finite limit for /? — > 00, or — > 0. 

For each /3, we have carried out simulations at a number of different lattice extents A^; 
examples are shown in Fig. El No significant volume dependence is observed for fi/N < 1, 
and we thus estimate the infinite-volume limit by fitting a constant to data in this range. 

Given the infinite-volume estimates, we extrapolate the data to the continuum limit, /? 
00. In Fig. H we show the functions - ^Tr [Pi2])a - [ci//3 + 02//?^ + cs//?^]} and 

- iTr [Pi2])a - [ci//3 + C2//32 + C3//33 + C4 In/?//?^]}. It is observed how even the 4-loop 
logarithmic divergence is visible in the data, as some upwards curvature for < 0.06. On 
the other hand, for < 0.01 the significance loss due to the subtractions grows rapidly and 
the error bars become quite large, so that these data points have little effect on the fit. 

The continuum extrapolation is carried out by fitting a function di + ^2 //? + to the 

infinite- volume extrapolated data for { ( 1 - ^ Tr [P12] ) a - [ci //? -|- C2 + 03//?^ -|- C4 In /J//?"^] } , 
in the range 0.01 < 1//5 < 0.10. We find that this functional form describes the data very well. 
The fitted values are di = 19.4.. .20.7, ^2 = 110.. .63, = 717.. .1101, with xVdof = 5.8/6, 
where the intervals indicate the projections of the 68% confidence level contour (i.e. the 
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Figure 2: The significance loss due to the subtractions of the various ultraviolet divergent 
contributions in the gluon condensate. Here again "plaq" = (1 — ^Tr [Pi2])a, and the symbols 
Cj in the curly brackets indicate which subtractions of Eq. ()2.12|) have been taken into account. 

surface where = Xmin + 3.53) onto the various axes, from one end of the elongated ellipse 
to the other. ^ We have also estimated the systematic errors from the effect of higher order 
terms in the fit ansatz, and found that they are of the same order as these intervals, which 
we thus consider as our combined error estimates. Returning back to Eq. (|2.12|) . we then 
obtain our final result, 

where we have inserted A'c = 3. 
4. Conclusions 

The purpose of this paper has been to study the expectation value of the elementary plaquette 
in pure SU(3) lattice gauge theory in three dimensions, as well as to outline how the MS 
scheme gluon condensate of the continuum theory can be extracted from it. To achieve 
this goal, we have carried out high precision numerical Monte Carlo simulations close to the 
continuum limit, corresponding to lattice spacings 0.05 < agf < 0.5, where is the gauge 
coupling. 

^If the three largest /3's are included in the fit, the parameters remain essentially the same, di = 19. 4. ..20. 8, 
d2 = 107.. .62, da = 733. ..1117, while x^/dof = 7.0/9 has decreased due to the large error bars at these /3's. 
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Figure 3: Finite-volume values for - ^Tr [Pi2])a - [ci//3 + 02//?^ + cs//?^ + C4 InP/P'^]}, 

as a function of the physical extent j3/N = G/g^L of the box. The solid symbols indicate the 
infinite- volume estimates, obtained by fitting a constant to data in the range fi/N < 1. 



When the leading perturbative terms, up to 4-loop level, are subtracted from the plaquette 
expectation value, and the result is divided by a finite quantity remains (the right-hand 

side of Eq. H2.12() . without C4) which can be taken as the definition of a renormalised gluon 
condensate in lattice regularization (in certain units). We have carried out the subtractions 
and the extrapolation ag^ 0, and shown that our data appear to be precise enough to 
determine the remainder with less than 5% errors, cf. Fig. and Eq. (|3.H) . 

To relate this number to the gluon condensate in some continuum scheme, say MS, a further 
perturbative 4-loop matching computation remains to be completed, fixing the constant C4 in 
Eqs. (|2.12j) . (|3.1|) . Our study should provide a strong incentive for finalising this challenging 
but feasible task, and there indeed is work in progress with this goal. The MS scheme 
conversion is also needed in order to apply our result in the context of finite temperature 
physics, particularly for determining the 0{g^T'^) contribution to the pressure of hot QCD, 
since the other parts of that computation have been formulated in the MS scheme |12| . 
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Figure 4: The infinite- volume extrapolated data, plotted as in Fig. |21 The effect of the 4-loop 
logarithmic divergence is to cause additional upwards "curvature" in the upper data set. The 
lower set includes all the subtractions, and should thus have a finite continuum limit. The 
continuum extrapolation (as described in the text) is indicated with the dashed line. The 
gray points have error bars so large that they are insignificant as far as the fit is concerned. 
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